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The formalism for neutrino flavor change induced by lepton family number violating interac¬ 
tions is developed for the three-neutrino case, and used to derive the corresponding flavor change 
probabilities in matter. Applied to the solar and atmospheric neutrino fluxes, it is argued that the 
observed anomalies, including the zenith dependence for the atmospheric case, could be due to such 
interactions. 
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I. INTRODUCTION 

The deficit in the observed solar neutrino fluxes (see for example reference Q), and the deviation of the atmospheric 
“ratio of ratios” from 1.0 [2j are two of the outstanding problems in today’s particle physics. These are usually 
interpreted as being due to neutrino oscillations, for which, however, there is as yet no unambiguous evidence. 

In this paper, the possibility that both the solar and atmospheric neutrino anomalies could be due to the existence of 
lepton number violating interactions, and in particular supersymmetric R-parity violating interactions, is investigated. 
Whereas various authors || |(J have recently suggested these interactions could lead to neutrino masses and mixings 
leading to the “standard” vacuum oscillation solution, here possible direct, matter-induced flavor change is studied. 


II. FLAVOR CHANGING INTERACTIONS AND THE MSW EFFECT 


The time dependency of the flavor composition of neutrinos propagating in matter is governed by the equation 
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where the term \/2 GpWe - N e is the electron density - corresponds to the usual v e — e _ charged-current interactions 
(this is the term at the origin of the classic MSW effect Jt[[| ), and the cu and bij coefficients represent flavor-diagonal 
(FD) and flavor-changing (FC) interactions, respectively. Note that we have assumed CP conservation here (i.e. the 
propagation matrix is symmetric). If such interactions exists, they would lead to additional electroweak potential 
energy for the propagating neutrino eigenstates, thereby opening the possibility for resonant enhancement of flavor 
change. In the Standard Model, electroweak neutral currents are FD interactions, while FC interactions do not exist 
in the leptonic sector. 

In supersymmetric theories with R-parity violation however, we find both new FD and FC interactions may exist. 
The two types of interactions that are relevant for neutrinos give rise to v — d-type quark and v — e scattering via 
d-type squark and selectron exchange, respectively. They correspond to the following two terms in the lepton number 
violating part of the superpotential: 


W = A,,,/.',/.'+ X[ jk L i L Q j L D k R . (2) 

Here i,j and k are the generation indices, Lr and Qr are the left-handed lepton and quark doublet superfields, and 
Er and Dr are the right-handed lepton and quark singlets. Note that the A couplings are antisymmetric in the first 
two indices. 

For example, for neutrino energies small compared to the squark mass ( rriq ), the term 
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corresponds to flavor-changing neutrino scattering on a d-quark where the &-squark acts as the mediator of the 
interaction. We can now rewrite the coefficients of the propagation matrix in equation (ij) in terms of the R-parity 
violating couplings. Taking only into account scattering on d quarks (matter is essentially first-generation) we find 
for example: 
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where N d is the down-quark density in matter. All the other matrix coefficients may be expressed likewise. 

Various authors have suggested @0 or established |yj that the solar neutrino problem could be explained by the 
presence of non-standard FD and FC interactions. The solutions have recently been analyzed again by P.I. Krastev 
and J.N.Bahcall J0] taking into account the newly available solar neutrino data. 

However, all these studies only considered the solar neutrino problem, and used a simplified two-family model 
(involving the first and third generations). In this work, a 3-family analysis is developed, in the attempt to explain 
both the solar and atmospheric neutrino anomalies. 

The procedure to derive the flavor change probabilities in the 3-family case is as follows. Let \u a > be the flavor 
eigenstates (a = e, /r or r) and \vi > the propagation eigenstates (*=1,2 or 3). We then have 


3 3 3 / 3 \ 
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where U a i is a unitary rotation and the are the propagation energy eigenlevels, eigenvalues of the propagation 
matrix given in equation (|I]). Since the elements of the propagation matrix, hence also its eigenvalues, depend on 
matter density it is advisable to replace the time t by the distance traveled x with the approximation t = x. The 
flavor changing amplitude is then given by 
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However, since the columns of the rotation matrix U which diagonalizes the propagation matrix are precisely the 
eigenvectors of this matrix,one finally obtains 


P (u a (x = 0)-* vp(x = L)) = 
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where V is the matrix whose columns are the eigenvectors of the propagation matrix. 


( 7 ) 


III. DATA 


Combining solar and atmospheric neutrino data yields 13 constraints: there are three different types of solar neutrino 
experiments, each sensitive to a different energy range, and the atmospheric “ratio of ratios” is measured in 5 zenith 
angle bins for two different energy ranges (sub- and multi-GeV). 

The three types of solar neutrino experiments (Gallium: Sage and Gallex, Chlorine: Homestake, and water 
Cerenkov: Super-Kamiokande) have different lower thresholds for the neutrino energy. This has allowed to show 
that there seems to be an energy-dependent suppression of the solar neutrino flux (see for example N. Hata and P. 
Langacker jl]]), which could be a strong hint for neutrino oscillations, were it not that neutrinos of different energies 
are prod uce d at different radii in the Sun and therefore “see” different matter densities on their way out. 

Table III summarizes the results from the various experiments (the results from the Gallium experiments have 
been combined into one), and gives the ratio of observed rates over rates predicted by the Bahcall-Pinsonneault 95 
Standard Solar Model |h|. The experimental results can be found in Refs. JmJ] (Homestake), Jh| (Sage), jT(| (Gallex) 
and |l7j (Super-Kamiokande). 
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For the atmospheric neutrino data, the results published recently by the Super-Kamiokande Collaboration |2| are 
being used. Fogli et al. jl 8 ) have graphically extracted the numerical values from the Super-Kamiokande plots, yielding 
the observed number of data events and the expected number of events both for data and Monte-Carlo, such that all 
the detection efficiencies can be taken into account. Table III gives the observed “ratio of ratios” (ratio of muon-like 
over electron-like events for data divided by the same ratio for Monte-carlo) for each of the five zenith bins, both for 
the sub- and multi-GeV samples. The first error is statistical (ll| and the second is systematic Q. 

We have extracted upper limits on the values of the coefficients Cu and bij from the present limits [^9|,|o| on R- 
parity violating couplings. For reasons of simplicity, only the couplings (corresponding to scattering by squark 
exchange) and scattering on d quarks are included, so that the limits are given by equations analogous to equation 
( 0 ). These limits have then been transformed to limits on 6 variables en, £ 12 , 613 , £ 22 , £23 and £33 defined by: 
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which give the relative strength of the interaction w.r.t. the weak interaction. 

Furthermore, we can subtract cu x 1 from the right-hand side of equation (|lj). This does not affect our result since 
it is equivalent to a redefinition of the zero of the neutrino electroweak potential energy. 


IV. ANALYSIS AND RESULTS 
A. Analysis Procedure 

To find the values of the couplings which would produce effects similar to those observed in the data, a y 2 -like 
function with 13 terms corresponding to the three types of solar neutrino experiments and the atmospheric data 
is constructed. The minimum of this y 2 is then searched for (using Minuit [jiTf ), as a function of the 4 variables 
£ 13 ) 623,622 and £ 33 . £12 has been excluded from the fit since the strong limits arising from the non-observation of 
pTi^eTi ]l9| make its contribution to any of the two anomalies negligible. 

The Solar model used is the Bahcall-Pinsonneault Standard Solar Model |l3), which “cuts” the Sun into layers, 
and for each of these layers, gives the average electron and neutron densities (the Sun is assumed neutral so that the 
proton density equals the electron density) as well as the fraction of neutrinos produced for each of the major neutrino 
sources (pp, 8 B , 13 iV, ls O , 7 Be and pep). Our program steps through each layer, adds the newly produced neutrinos 
and changes the flavor of a fraction of the entering and produced neutrinos in the manner determined by the flavor 
change probability Eq. [i], which is a function of the variables and the densities. 

The Earth is approximated by 5 concentric spheres, each of uniform density and composition as has been done 
by Giunti et al. ^2|. To find the average distance that neutrinos travel through each layer as a function of their 
momentum and zenith angle, a small Monte-Carlo simulation has been developed. The simulation assumes a flat 
zenith angle distribution in each bin, smears it with the average angular correlation between the charged lepton and 
the neutrino (55 and 20 degrees for sub-GeV and multi-GeV neutrinos, respectively [(|), and then tracks the neutrinos 
back through the Earth, computing the distance traveled through each layer. The values found are given in Table |. 

The algebraic expressions for the eigenvectors and eigenvalues of the propagation matrix as a function of its coeffi¬ 
cients are fairly complex. We have used Mathematica P3[ to generate the expression in computer-readable form and 
systematically checked that the orthogonality conditions are satisfied. This is done to verify that solutions are not 
produced artificially due to machine accuracy. 


B. Results 


Table |I| gives the expected value for 
value found is 15.6, corresponding to 


The values of the parameters resulting in the lowest y 2 are given in Table [n] 
each of the experimental results for the best fit parameters. The minimum 
a probability of 7.6%. 

The variables are somewhat correlated (as expected) such that changes in their values generally affect the width of 
the “good” region for the others. Figure [l] shows the region in parameter space allowed at 95 % and 99 % C.L. for 
each combination of off-diagonal versus diagonal parameter. 

Since the value found for £23 has to be significantly larger than what is allowed by the present non-observation of 
the decay r —> p°p, the new interactions must either only apply to neutrinos, or suffer significant SU(2)l violation. 
It has recently been argued |24[|, however, that such violation could be in excess of present experimental bounds. 
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V. CONCLUSIONS 


We have searched for, and found coupling values of a lepton family number violating interaction which would explain 
both the solar and atmospheric neutrino anomalies through an MSW-like mechanism in the Sun and in the Earth. 
The fit to the data is not as good as for a purely neutrino oscillation solution, and the interaction must either only 
apply to neutrinos or undergo significant SU{2)l violation. It is however clear that such an interaction could play a 
major role in explaining the anomalies. 

It should be noted that in order to have significant flavor change, neutrinos have to traverse large amounts of 
matter. This results in a large predicted difference between the atmospheric ratio of ratios for the first two and last 
three bins of the Super-Kamiokande multi-GeV azimuthal dependence histogram. Increased statistics may make this 
effect visible. 
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FIG. 1. Allowed regions in parameter space for a) £13 versus £ 22 , b) £23 versus £ 22 , c) £13 versus £ 33 , and d) £23 versus £ 33 . 
Both 95 and 99 % C.L. regions are indicated. 


TABLE I. Average distance (in km) traveled by neutrinos through each Earth layer (l=core) for each zenith and momentum 
bin. 


Sample 

Bin No. 

cos 6 

Layer 1 

Layer 2 

Layer 3 

Layer 4 

Layer 5 

Sub-GeV 

1 

[- 1 . 0 , - 0 . 6 ] 

245.0 

1780.8 

3663.3 

630.4 

1131.5 


2 

[- 0 . 6 , — 0 . 2 ] 

149.4 

1130.1 

2623.4 

505.4 

983.7 


3 

[- 0 . 2 , 0 . 2 ] 

81.8 

649.9 

1731.3 

373.3 

791.1 


4 

[ 0 . 2 , 0 . 6 ] 

37.6 

318.5 

1000.0 

243.8 

565.7 


5 

[ 0 . 6 , 1 . 0 ] 

10.9 

100.9 

396.4 

112.1 

291.6 

Multi-GeV 

1 

[- 1 . 0 , - 0 . 6 ] 

247.5 

2148.3 

5157.8 

811.3 

1237.8 


2 

[- 0 . 6 , — 0 . 2 ] 

8.2 

236.4 

2398.9 

775.0 

1671.9 


3 

[- 0 . 2 , 0 . 2 ] 

0.1 

12.0 

466.8 

286.5 

1029.4 


4 

[ 0 . 2 , 0 . 6 ] 

0.0 

0.2 

34.1 

36.1 

237.7 


5 

[ 0 . 6 , 1 . 0 ] 

0.0 

0.0 

0.4 

0.8 

11.7 
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TABLE II. Fit variables, maximum allowed absolute value derived from the present limits on R-parity violating couplings 
(assuming SU{2)l symmetry), and lowest results. The corresponding y 2 value is also given. 


Variable 

Maximum Absolute Value 

Best Fit Value 

ei3 

0.421.ICE 1 

0.45.10“ 2 

623 

0.708.10~ 2 

0.19 

622 

0.164 

0.05 

633 

0.460 

0.414 

X 2 


15.6 


TABLE III. Fit results compared with the data. The first column indicates each of the 13 constraints and the second one 

gives the observed results: 

for the solar data the fraction of observed rate over 

the rate predicted by the Bahcall- Pinsonneault 

Standard Solar Model 95, and for the atmospheric data the “ratio of ratios” 
column shows the predicted results for the model described in this paper. 

(see text) for each zenith angle bin. The third 

Solar Neutrinos 
Experiment Type 

Observed Flux Fraction 

Best Fit Flux Fraction 

Homestake 

0.273 ±0.021(ea:p) t° 0 °l (theo) 

0.368 

Gallium 

0.507 ± 0.049(ea:p) ± 0.0‘28{theo 

) 0.424 

Sup er- Kamiokande 

0-366 ±°;”I iexp) t° n °^(theo) 

0.372 

Atmospheric Neutrinos 
Sub-GeV ( cos6 range) 

Observed Ratio of Ratios 

Best Fit Ratio of Ratios 

[-1.0,-0.6] 

0.41 ± 0.07(stat) ± 0.08(sj/st) 

0.50 

[—0.6, —0.2] 

0.63 ± 0.07 (stat) ± 0.08(s2/st) 

0.50 

[—0.2, 0.2] 

0.63 ± 0.07 (stat) ± 0.08(sj/st) 

0.58 

[0.2, 0.6] 

0.82 ± 0.07 (stat) ± 0.08(sj/st) 

0.79 

[0.6,1.0] 

0.71 ± 0.07 (stat) ± 0.08(sj/st) 

0.96 

Multi-GeV ( cosd range) 

Observed Ratio of Ratios 

Best Fit Ratio of Ratios 

[-1.0,-0.6] 

0.42 ± 0.U{stat) ± 0.12 (syst) 

0.51 

[-0.6, —0.2] 

0.52 ± 0.13(stat) ± 0.12 {syst) 

0.54 

[—0.2, 0.2] 

0.70 ± 0.12(stat) ± 0.12(sj/st) 

0.90 

[0.2, 0.6] 

0.72 ± 0.12(stot) ± 0.12(syst) 

1.00 

[0.6,1.0] 

0.87 ± 0.16(stat) ± 0.12(syst) 

1.00 
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